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Abstract. We present observations of 11 high redshift quasars (3.9 S z :£ 5.0) observed with low spectral resolution in the 
restframe ultraviolet using FORS 1 at the VLT UT 1. The emission-line fluxes of strong permitted and intercombination ultravi¬ 
olet emission lines are measured to estimate the chemical composition of the line emitting gas. Comparisons to photoionization 
calculations indicate gas metallicities in the broad emission line region in the range of solar to several times solar. The average 
of the mean metallicity of each high-z quasar in this sample is Z/Zq = 4.3 ± 0.3. Assuming a chemical evolution time scale 
of Tjvd; - 0.5 - 0.8 Gyrs, we derive a redshift of z/ ^ 6 to 8 for the onset of the first major star formation episode (Ho = 65 
km s“* Mpc^', Qm = 0.3, = 0.7), corresponding to an age of the universe of several 10* yrs at this epoch. We note that this 

epoch is also supposed to be the era of re-ionization of the universe. 
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1. Introduction 

Quasars are among the most luminous objects known in the 
universe. Due to their high luminosity they are excellent tools 
to probe their galactic environment up to the highest redshifts. 
There is growing evidence that quasar activity and the forma¬ 
tion of their host galaxies, in particular of massive spheroidal 
systems, are closely related. The detection of dark massive ob¬ 
jects (DMOs) in the center of nearly every galaxy with a sig¬ 
nificant spheroidal subsystem provides further strong support 
for the relationship between galaxy formation and quasar ac¬ 
tivity (Kormendy & Richstone 1995; Magorrian et al. 1998; 
Kormendy & Gebhardt 2001). It has been shown that the mass 
of the DMOs, most likely supermassive black holes, is closely 
correlated to the bulge mass of the host galaxy (Gebhardt et 
al. 2000; Kaspi et al. 2000; Merritt & Ferrarese 2001; Tremaine 
et al. 2002). The evolution of the star formation rate indicates 
that it was more than one order of magnitude larger at epochs 
z ^ 1 than in the local universe (Lilly et al. 1996; Connolly 
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et al. 1997; Tresse & Maddox 1998; Steidel et al. 1999). Strong 
evidence for the relationship of quasar activity to galaxy forma¬ 
tion accompanied by intense star formation is provided by the 
detection of large amounts of dust (~ 10* Mq) and molecular 
gas (~ 10'° Mq) in high redshift quasars (Andreani, La Franca, 
& Cristiani 1993; Isaak et al. 1994; Omont et al. 1996,2001; 
Carilli et al. 2001). Several galaxy evolutionary models have 
been suggested which show that galactic spheroids can eas¬ 
ily reach solar or supersolar gas-phase metallicities on time 
scales shorter than ~ 1 Gyr (Arimoto & Yoshii 1987; Hamann 
& Ferland 1993; Gnedin & Ostriker 1997; Friaga & Terlevich 
1998; Cen & Ostriker 1999; Salucci et al. 1999; Kauffmann & 
Haehnelt 2000; Granato et al. 2001; Romano et al. 2002). As a 
result, quasars at high redshift are valuable probes for dating 
the first star formation in the early universe. In particular, z ^ 4 
quasars probe a cosmic era when the universe had an age of 
less than ~ 10% of its current age (assuming Ho - 65 kms^' 
Mpc-', Q.M = 0.3, Qa = 0.7). 

The prominent emission-line spectrum of quasars con¬ 
tains valuable information to estimate the gas metallicity 
at early epochs due to star formation (for a review, see 
Hamann & Ferland 1999). Early studies to estimate the 
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abundances in broad emission line region (BELR) gas were 
based on several generally weak intercombination lines like 
Niv]/I1486, Oiii]/ll663, Niii]dl750, and Ciii]dl909 (Shields 
1976; Davidson 1977; Baldwin & Netzer 1978; Osmer 1980; 
Gaskell, Shields, & Wampler 1981; Uomoto 1984) and indi¬ 
cated already higher than solar metallicity for the BELR gas. 
Recent studies of the emission line and intrinsic absorption line 
properties of (z - 3) quasars provide evidence for enhanced 
metallicities up to an order of magnitude above solar (Hamann 
& Eerland 1992,1993; Petitjean et al. 1994; Eerland et al. 1996; 
Hamann 1997; Pettini 1999; Dietrich et al. 1999; Dietrich & 
Wilhelm-Erkens 2000; Hamann et al. 2002; Warner et al. 2002; 
Dietrich et al. 2002 in prep.). These high metallicities require a 
preceding intense star formation phase. 

To estimate the chemical composition of the gas in quasar 
BELRs, nitrogen as a secondary element is of particular in¬ 
terest. Providing that the secondary nitrogen production, i.e., 
the synthesis of nitrogen from existing carbon and oxygen via 
CNO burning in stars (Tinsley 1979,1980; Wheeler, Sneden, 
& Truran 1989), is the dominant source for nitrogen, we can 
expect the relation N/0 ~ OjH ~ Z, i.e., NjH ~ 7}. This 
scaling of NjH with metallicity has been confirmed for many 
Hii-regions (Shields 1976; Pagel & Edmunds 1981; van Zee et 
al. 1998; Izotov & Thuan 1999). 

As suggested by Shields (1976), and later developed by 
Hamann & Eerland (1992,1993) and Eerland et al. (1996), 
emission line ratios involving Nv/I1240 are especially valu¬ 
able. Generally, Nv/I1240 is stronger than expected in 
the spectra of high redshift quasars compared to pre¬ 
dictions of standard photoionization models assuming so¬ 
lar metallicity. Hamann & Eerland (1992,1993) showed 
that N vT 1240/C ivT 1549 and N vT 1240/He ii/l1640 are use¬ 
ful metallicity indicators. Recently, Hamann et al. (2002) 
presented results of a detailed study of the effects of metal¬ 
licity and the spectral shape of the photoionizing contin¬ 
uum on emission line ratios. They revised the metallicity de¬ 
pendence of line ratios involving intercombination lines, re¬ 
combination lines, and collisionally excited lines. They sug¬ 
gest that the most robust indicators of the gas chemical 
composition are the line ratios Niii]/I1750/Oiii]/I1663 and 
N vT1240/(O vi/11034-hC ivT1549). 

We observed a small sample of high redshift quasars (z ^ 4) 
with the Very Large Telescope (VLT) to extend prior studies to 
higher redshift and hence to earlier epochs in the cosmic evolu¬ 
tion, approaching an age of the universe of ~ 1 Gyr. This sam¬ 
ple is supplemented with the observation of SDSS 0338H-0021 
(z - 5.0; Dietrich et al. 1999; Ean et al. 1999). 

In section 2 we describe the observations and the data anal¬ 
ysis. In section 3 we present the results of the analysis of the 
emission line spectra. We estimate the elemental abundance 
of the line emitting gas based on the line ratios of several 
diagnostical emission lines (Hamann et al. 2002). The over¬ 
all mean metallicity of the high redshift quasars amounts to 
Z/Zo = 4.3 + 0.3. These results are discussed and are compared 
with recent studies in section 4. The chemical composition of 
the BELR gas provides further evidence that the first major star 
formation epoch started at a redshift of z/ - 6 to 8, correspond¬ 
ing to an age of the universe of several 10* yrs. This result is in 


Table 1. The high-z quasar sample 


quasar 

m-r 

z 

Q 0046-293 

19.38 

4.01 

Q 0101-304 

20.06 

4.07 

SDSS 0338+0021 

21.7 

5.00 

PKS 1251-407 

19.90 

4.46 

APM 1335-0417 

19.40 

4.38 

BRl 1500+0824 

18.84 

3.95 

BRI 1557+0313 

19.80 

3.89 

Q 2133-4311 

20.85 

4.18 

Q 2133-4625 

20.96 

4.18 

Q 2134-4521 

20.15 

4.37 

PC 2331+0216 

19.98 

4.09 


good agreement with recent model predictions relating quasar 
activity with the formation of early type galaxies. We assumed 
Ho = 65 km s ' Mpc ', = 0.3, and Da = 0.7. 


2. Observations and Data Analysis 

We observed the restframe wavelength region ~ 800 - 2000 A 
for 11 quasars with z ^ 4 (Table 1), which contains the valu¬ 
able diagnostic ultraviolet emission lines (e.g., O vid 1034, Lya, 
Nvdl240, Niv]dl486, Civdl549, Heiidl640, Oiii]T1663, 
Niii]/I1750, and C iii]/ll909). The quasars were observed on 
Dec. 26, 1998, July 16 - 20, and August 14 - 15, 1999. 
The total exposure times are ranging from 15 to 60 minutes, 
typically 30 minutes. All observations were carried out us¬ 
ing EORS 1 (focal reducer and low-dispersion spectrograph; 
Mohler et al. 1995) at the Cassegrain focus of the VLT UT 1 
Antu. The observations were performed under non-photometric 
conditions. The seeing was ~ 1" on December, and varied 
between O'.'8 and 2'.'6 during the observations in July and be¬ 
tween T.'O and T.'9 in August. A Tektronix CCD detector with 
2048x2048 pixels (pixel size lApm'KlApm) and grating G150I 
(230 A/mm) were used in connection with a blocking filter 
GG435 to suppress contamination of the spectra by the second 
order for T ^ 8700 A. With this setting an observed wavelength 
range of approximately 4200 - 10500 A was achieved. The 
quasar spectra were recorded in the multi-object spectroscopy 
(MOS) mode, in the 10'* (1" x 20") or 11'* (1" x 22") of the 
19 slitlets. The position angle of the slit was set perpendicular 
to the horizon to minimize light losses caused by differential 
refraction. With the exception of ARM 1335-0417, the quasars 
were observed close to the meridian (airmass less than ~ 1.1). 
The other slits were used to observe spectra of faint objects in 
the close environment of the quasars which appear slightly ex¬ 
tended on the Palomar Observatory Sky Survey (POSS). These 
serendipity data will be presented and discussed in an upcom¬ 
ing paper. The standard stars EG21, LTT 1788, and LTT7379 
(Hamuy et al. 1992) were observed for relative flux calibration 
each night, with exposure times of typically 60 seconds. 
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The quasar and standard star spectra were processed using 
standard MIDAS[| software. Cosmic-ray events were removed 
manually by comparing multiple exposures for each object. 
The night sky component of the 2D-spectra was subtracted 
by fitting third order Legendre polynomials perpendicular to 
the dispersion, along each spatial row of the spectra using ar¬ 
eas on both sides of the object spectrum which were not con¬ 
taminated by the quasar or other objects. The 1 D-spectra were 
extracted using an optimal extraction algorithm (Horne 1986). 
The width of the spatial extraction windows was adjusted to 
match the seeing recorded during the observation. The helium- 
argon wavelength calibration frames, taken for each quasar 
MOS-setting, yield a dispersion of 5.4A/pxl with an internal 
error of Ad ^ 0.4 A. The strong night sky emission lines [O i] 
45577, [O i] 46300, and [O i] 46364 indicate an absolute un¬ 
certainty of A4 ^ 1.1 A (Av ^ 50kms"'). The FWHM spectral 
resolution measured from these lines is A4 ^ 25 A. 

We corrected each quasar spectrum for the atmospheric c- 
band, b-band, and A-band absorption, caused by O 2 , and the 
atmospheric water vapor absorption bands (447140 - 7340 A, 
8140 - 8350 A, 9250 - 9600 A), using the standard stars ob¬ 
served during the same nights. The spectra were corrected for 
atmospheric extinction applying the standard curve of La Silla 
(Schwarz & Melnick 1993) and for interstellar extinction using 
the Eb-v values of Burstein & Heiles (1982) and the extinction 
curve of Savage & Mathis (1979). 

The sensitivity functions, based on individual standard 
stars, differ by less than ~ 4 % from the mean sensitivity func¬ 
tion for all nights. 

3. Results 

The flux calibrated quasar spectra are shown in Figure 1, 
with the observed flux in units of 10“'^ ergs^' cm“^ A“*. The 
strongest lines in the spectra are Lya41216 and C iv41549. The 
broad and moderately strong 41400 A feature which consists 
of the Siiv41397,1403 and the Oiv]41402 multiplet is quite 
prominent in these quasar spectra. This feature will be ref- 
ered as Siiv41400 in the following. In addition to these emis¬ 
sion lines, several important diagnostic lines like Ovi41034, 
Heii41640, Ora]41663, Niii]41750, and Cra]41909 are visi¬ 
ble and marked in Figure 1. 

The observed quasar spectra were transformed to their rest- 
frame using the redshifts given in Table 1. To determine the 
redshift we fit a Gaussian profile to the upper part of the the 
Civ41549 emission line (/^ > 50% of the peak intensity). 

We employed a multicomponent fit to the quasar spectra to 
determine the power-law continuum, Fy oc v“, the contribution 
of Fe II and Fe iii emission, and the weak contribution of the 
Balmer continuum emission (see Dietrich et al. 2002a,b). The 
power-law continuum was fitted using small spectral regions, 
each 10 to 20 A wide, centered at 4 ^ 1290 A, 1340 A, 1450 A, 
1700 A, 1830 A, and 1960 A which are free of detectable emis¬ 
sion lines. The spectral indices a are given for each quasar in 
Table 2. The properly scaled Fe emission was subtracted using 

* Munich Image Data Analysis System, trade-mark of the European 
Southern Observatory 
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Fig. 2. In the top panel an example of the Lya, Nv41240 
emission line profile complex reconstruction is shown for 
PKS 1251-407. The individual profiles are displayed as solid 
lines while the broad and narrow components are shown as 
dotted lines. The sum of the profile fits of the individual 
lines is plotted as short dashed line. At the bottom of the 
figure we plotted the difference between the individual pro¬ 
file fits and the power-law continuum corrected quasar spec¬ 
trum. In the bottom panel, the reconstruction of the C iv41549, 
Heii41640, Oiii]41663 emission line profile complex is dis¬ 
played for PKS 1251-407. 


an emission template which accounts for both Fe ii and Fe iii 
emission (Vestergaard & Wilkes 2001). This improves espe¬ 
cially the measurement of N iii]41750 and reduces the flux of 
the 41600 A feature (Laor et al. 1994). 

We used the Civ41549 emission line profile as a template 
to measure the other emission line fluxes. The Civ41549 line 
profile was fitted with a broad and a narrow Gaussian compo¬ 
nent. We fixed the spectral width of the broad and narrow com¬ 
ponent in velocity space and allowed the strengths to vary in¬ 
dependently. Furthermore, shifts in velocity space of the broad 
and narrow component were restricted to a range of less than a 
few 100 kms^' with respect to C iv41549. Using this Civ41549 
emission template is particularly important for the Nv41240 
and Heii41640 profiles which are blended with other emis¬ 
sion lines. To measure C iii]41909 we used the C iv41549 tem¬ 
plate to fit Ciii]41909, Siiii]41892, and Aliii41857 simulta¬ 
neously. This template-fitting approach can be justified since 
Civ41549, Nv41240, and Heii41640 are all high ionization 
lines (HIL). Figure 2 shows typical examples of the deblend¬ 
ing of the Lya41216 - Nv41240 and Niv]41486 - Civ41549 
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- He 11/11640 - Oiii]/ll663 emission line complexes. The re¬ 
sulting line flux measurements for the quasars are given in 
Table 2. The uncertainties were estimated from the multi- 
component line fit using the scaled Civ/11549 line profile to 
obtain a minimum fit. The errors are of the order of ~ 
10 % for stronger lines like LyaAl2l6, N V/11240, Si iV/11400, 
C iV/11549, and C iii]/11909, and ~ 20 % or more for the weaker 
lines. The measurement of the N iv]/ll486 emission line flux is 
severely affected by the blue wing of the broad C iv-component. 
Particularly, for quasars with broad emission line profiles the 
N iv]/ll486 line tends to show a low contrast to the outer part 
of the C IV/11549 line profile which can be represented by the 
broad C iv-component only (see Dietrich & Hamann 2002 in 
prep., for further discussion). 

The O VI/11034 emission line flux given in Table 2 has 
been corrected for Lya forest absorption. For this correction 
we assumed an intrinsic continuum slope of a = -1.76 for 
A < 1200 A (Telfer et al. 2002). The correction factor follows 
by simply assuming the same fraction of the continuum and 
O VI/11034 emission line flux were absorbed by the Lya forest. 
The correction factors ranged from 1.33 to 2.67, with an aver¬ 
age of 1.71 +0.42. 

Hamann et al. (2002) presented new results of the depen¬ 
dence of emission line ratios from metallicities and the shape 
of the ionizing input continuum. They calculated detailed pho¬ 
toionization models for a wide range of density, n^, and con¬ 
tinuum strength, <pH, using CLOUDY (Ferland et al. 1998). To 
study the influence of the chemical composition and of the 
continuum shape, the metallicity was varied from Z/Zq = 0.2 
to 10 with three different input continua — a broken power- 
law continuum with a UV bump (Mathews & Ferland 1987), 
a single power-law continuum with a - -1, and a segmented 
power-law based on Zheng et al. (1997) and Laor et al. (1997). 
For more details of the model calculations see Hamann et 
al. (2002). We used the results of these calculations to derive 
the chemical composition of the line emitting gas. The most 
reliable line ratios are those of N iii]/0 iii] and N v/(0 vi+C iv). 
In general, the derived metallicities depend only weakly on 
the shape of the ionizing continuum. The strongest influence 
of the input continuum shape is found for N v/He ii, where we 
find a factor of 2 difference depending on the continuum. Our 
estimates of the metallicity for the high redshift quasars un¬ 
der study (Table 3) are based on model results using the seg¬ 
mented power-law continuum, mentioned above. The metallic¬ 
ities based on the other spectral shapes deviate ^ 25 % from the 
values derived for a broken power-law continuum. The largest 
deviations occur for the N v/He ii line ratio. This line ratio is 
quite sensitive to the temperature of the gas and on the ioniz¬ 
ing continuum shape. However, Hamann & Ferland (1993) and 
Ferland et al. (1996) showed that this line ratio provides a firm 
lower limit on N/He when adopting a BELR ionization param¬ 
eter that maximizes the predicted N v/He ii line ratio, e.g. a hard 
power-law continuum (a - -1.0). 

Ratios of the line fluxes in Table 2 were transformed to 
metallicity estimates of the gas using Figure 5 in Hamann et 
al. (2002). The chemical composition of the gas as provided by 
each of the emission line ratios is plotted as a function of red¬ 
shift in Figure 3. The most suited line ratios, Niii]/Oiii] and 


N v/(0 vi+C iv), as well as N v/O vi, N v/C iv, yield consistent 
metallicity estimates. The errors of the metallicities, given in 
Table 3, were derived from the errors of the line flux measure¬ 
ments. These uncertainties were used to calculate the uncer¬ 
tainty of the line ratios which yields the range of metallicities 
for a given line ratio, using the results presented in Hamann et 
al. (2002). 

In particular, gas metallicities based on N iii]/0 iii] and 
N v/(0 VI+C iv) agree within ~ 33 %. The ratio N iii]/C iii] tends 
to indicate lower metallicities than the other line ratios. This 
may be due to the fact that the ratio depends more on the tem¬ 
perature of the gas than N iii]/0 iii]. Furthermore, the critical 
densities of N iii] and C iii] differ by more than a factor of 2. 
Hence, this line ratio is not as robust as N iii]/0 iii], since uncer¬ 
tainties are introduced if the emission is received from spatially 
different parts of the line emitting region with different physi¬ 
cal conditions (Hamann et al. 2002). Overall we conclude that 
the line ratios involving N iii] and N v provide quite consistent 
estimates of the gas metallicity for quasars. 

For a few quasars we could also compare the metallici¬ 
ties based on N iv]/0 iii] with the results using ratios including 
N III] and N v, respectively. The metallicities we obtained using 
N iv]/0 III] are in agreement with the chemical composition of 
the gas based on N iii] and N v line ratios to within ~ 50%. 
However, the metallicities obtained by analyzing N iv]/C iv are 
lower than the other ratios (Table 3, Figure 3). The tendency 
to significantly lower metallicity estimates based on N iv]/C iv 
was also noted by Shemmer & Netzer (2002). This is presum¬ 
ably due mostly to the difficulty in measuring the line flux of 
the weak Niv]/11486 emission line. In particular, for quasars 
with broad emission line profiles the Niv]/H486 line is lo¬ 
cated in the outer wing of the Civ/11549 profile. With a typ¬ 
ical strength of about ~ 5 % of the Civ/11549 line flux for 
solar metallicities, this line can be well hidden in the outer 
blue wing of Civ. The physical reason for the discrepancy is 
not understood. However, it is important to keep in mind that 
this emission line ratio compares an intercombination line to a 
strong permitted line which may originate under different phys¬ 
ical conditions. 

For each quasar, we used all of the available metallicity es¬ 
timates based on the individual emission line ratios to calculate 
the mean metallicity. With the exceptions of SDSS 0338+0021 
and Q 2133-4625, 4 to 8 individual estimates were averaged. 
The resulting metallicities are typically several times solar 
(Table 3). The mean metallicity for each quasar, given as the 
average of the individual estimates, is shown as a function of 
redshift in Figure 4. The overall average gas chemical com¬ 
position for the 11 high redshift quasars, using these mean 
metallicity estimates, is Z/Z© = 4.3 + 0.3. [| We also calcu¬ 
lated the mean metallicity given by each line ratio (Table 3, 
last line). The most robust metallicity indicators, Niii]/Oiii] 


^ Recent studies on the photospheric solar abundance indicate car¬ 
bon and oxygen are about 30% lower than the values given by 
Grevesse & Sauval (1998), while the nitrogen abundance remains 
nearly unchanged within the errors (Holweger 2001; Allende Prieto 
et al. 2001). Hence, this might result in a reduced estimate of the su¬ 
persolar abundances for the high redshift quasars by 30 %. 
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Table 2. Broad emission-line flux measurements for the observed high redshift quasars. The integrated restframe emission line 
fluxes are given in units of 10“^^ergs“^ cm“^ and L^(1450) in ergs^^ In addition, the spectral index a of the power-law 
continuum fit is given (Fy ~ v“). 





T/,„[10-‘5 

erg s ' cm ^] 




Q0046-293 

Q0101-304 

SDSS0338±0021 

PKS 1251-407 

APM 1335-0417 

BRI 1500±0824 

log LA 1450) 

43.60 

43.59 

43.45 

43.32 

43.56 

43.49 

Ovid 1034 


289.0 ± 26.2 


82.7 ±9.8 


62.3 ± 15.6 

LyQ'1216 

505 ± 88 

850 ± 264 

290 ± 75 

282.6 ± 42.4 

480 ± 96 

360 ±46.1 

Nvdl240 

90± 18 

210 ±53 

83.1 ±20.8 

65.6 ± 8.8 

85± 17 

50.0 ±7.5 

Siiidl263 

25 ±2 

35 ±9 

6.0 ± 1.5 

18 ±3.0 

12.0 ±2.5 

13.0 ±2.0 

001305 

22.4 ± 0.6 

28.4 ±7.1 

13.0 ±2.5 

9.9 ± 1.0 

10.5 ± 1.5 

1L0± 1.8 

C11/11335 

15.2 ± 0.7 

16.2 ±4.1 

15.0 ±3.1 

2.0 ± 0.9 

8 ± 1.4 

8.0 ± 1.7 

SiiV/11402 

79.4 ± 8.5 

109.8 ± 22.0 

60.0 ± 8 

25.7 ± 1.3 

53.5 ± 9.8 

45.3 ±5.8 

Niv]/11486 


13.7 ± 1.8 


5.3 ± 1.6 

12.0 ±2.5 

6.3 ± 1.4 

CiV/11549 

163.3 ± 28.6 

276.1 ±85.7 

170 ± 43 

120.0 ± 9.4 

143.5 ± 29.6 

101.8 ± 13.0 

He 11/11640 

25.5 ± 2.2 

45 ±9 

40 ±8 

22.0 ±3.0 

27.9 ± 5.6 

17.0 ± 3.4 

O III] /H 663 

17.9 ± 1.9 

33 ±7 


7.0 ± 2.0 

17.5 ±3.5 

9.5 ± 2.6 

Niii]/H750 

18.0 ±2.7 

40 ±8 




12.5 ± 2.5 

A1III/H857 

27.0 ±4.1 





10.0 ± 2.0 

Siiii]/H892 

10.0 ± 2.0 

20 ±5 




21.8 ±2.9 

Ciii]/H909 

65.0 ± 9.9 

170 ± 35 




65.5 ± 12.8 

a 

0.09 ± 0.06 

-0.24 ±0.14 

-0.01 ± 0.20 

-0.11 ±0.17 

-0.45 ± 0.07 

-0.57 ± 0.06 


BRI1557±0313 

Q2133-4311 

Q2133-4625 

Q 2134-4521 

PC 2331±0216 


logLA1450) 

43.21 

42.82 

42.78 

43.15 

43.43 


OVI/H034 

83.5 ± 14.2 

34.8 ±6.3 

7.38 ± 1.40 

58.1 ± 12.2 

130.6 ± 25.7 


LyQ'1216 

275 ± 20.6 

166 ± 25 

18.9 ±5.6 

225 ± 45 

415 ± 62 


NV/H240 

45.5 ± 6.8 

27.7 ± 5.5 

3.4 ± 0.5 

34.2 ± 6.8 

89 ± 13.4 


Siii/H263 

2.0 ± 0.6 

4.7 ± 0.9 

1.3 ±0.3 

7.0 ± 1.4 

37 ± 6.5 


0 01305 

17.5 ± 2.6 

2.5 ± 0.4 


9.0 ± 1.8 

20.5 ±3.1 


C11/11335 





7.5 ± 1.2 


Siiv/H402 

10.3 ± 1.4 

5.5 ±0.8 

2.3 ± 0.6 

35.4 ±4.8 

42.0 ± 6.5 


Niv]/11486 

5.7 ± 1.2 

3.6 ± 1.0 


8.1 ±2.1 

3.4 ± 1.2 


CIV/H549 

128.8 ±9.7 

62.9 ± 9.5 

8.14 ±2.49 

85.9 ±21.5 

88.2 ±9.2 


He 11/11640 

30.0 ± 4.5 

13.4 ±4.0 


11.4 ±3.4 

13.5 ±2.0 


O III] /H 663 

14.0 ±2.1 

5.8 ± 1.8 


8.0 ± 2.4 

9.0 ± 1.4 


N III] /H 750 

14.5 ± 2.2 




8.3 ± 1.0 


Al 111/11857 

8.0 ± 1.5 






Siiii]/H892 

10.0 ± 1.7 




14.3 ± 2.2 


Ciii]/11909 

69.5 ± 10.5 




42.7 ± 6.4 


a 

-0.12 ±0.05 

-0.11 ±0.20 

0.71 ±0.47 

-0.15 ±0.23 

-0.74 ± 0.04 



and N v/(0 viH-C iv), infer quite consistent metallicities for this 
small sample of high-z quasars. Although the less suited line 
ratios, N iii]/C iii] and N iv]/C iv, indicate a lower metallicity, it 
is still at least solar. 

The emission line ratios indicate that the gas chemical com¬ 
position of the BELR at redshifts z ^ 4 is several times solar 
and at least of solar metallicity. In an ongoing study we are 
investigating the reason for the lower metallicities inferred us¬ 
ing N iv]/C IV and N iii] /C iii] compared to the robust line ratios 
N iii]/0 III] and N v/(0 vi-hC iv) to obtain solid measurements of 
the metallicity. 


4. Discussion 

The supersolar metallicities we derived for the emission line 
gas in quasars at z ^ 4 provide valuable information about 
the preceding star formation epoch. This early star formation 
epoch may well correspond to the beginning of major star for¬ 
mation in the host galaxies. In the context of one-zone chem¬ 
ical evolution models (e.g., Hamann & Ferland 1992,1993; 
Padovani & Matteucci 1993; Matteucci & Padovani 1993), as 
well as in more recent multi-zone models (Gnedin & Ostriker 
1997; Friaga & Terlevich 1998; Granato et al. 2001; Romano et 
al. 2002) which take into account the feedback of the dynami¬ 
cal and chemical evolution of the forming galaxies, super-solar 
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Table 3. Estimated relative abundance of the line emitting gas given in units of solar metallicity Zq. The estimates based on the 
results using a segmented power-law continuum fit for the photoionization models (Hamann et al. 2002). 


quasar 



Z/Zo 







Niii]/Oiii] 

Niii]/Ciii] 

Niv]/Civ 

Niv]/Oiii] 

N v/He II 

N v/C IV 

N v/O VI 

Nv/(Ovi±Civ) 

mean 

Q 0046-293 

3 

•^•^-0.6 

1 O+0.3 



6 4+2-^ 

-2.0 

5.4)];^ 



4.3 ±0.7 

Q 0101-304 

4 3+1-2 

1 

^•'^-0.3 

1 4+0.4 

^'^-0.4 

3-2-1:6 

10.3)1] 

7.9)];I 

4 3+1-^ 

^•*^-1.3 

6.1)|« 

4.9 ±0.9 

SDSS 0338-r0021 





-1.9 

4-7)?:9 



4.0 ± 1.4 

PKS 1251-407 



1 9-1-O.4 
^'^-0.4 

6.9!|^ 

4 7-^2.0 
-0.7 

5 4+*'^ 

-1.0 

4 7+1.8 
^•'-1.0 

4 Q+1.6 
^•^-1.1 

4.6 ±0.6 

ARM 1335-0417 



2 - 6 t^ 

6.3!|? 

5 0+^‘ 

5-8-?:8 



4.9 ±0.9 

BRI 1500-1-0824 

4 7+1-8 



6.U|> 

4 n+t.i 

-1.3 

4.7)1;] 

4 9+2.9 

^•^-1.6 

4 6+2-2 

^•^-1.6 

4.1 ±0.6 

BRI 1557-r0313 

-0.7 

\ 4+0.2 
-0.2 

1 9-I-O.4 
^'^-0.4 

3 0+1-^ 

9 1 4-0.6 
^'•^-0.4 

O 94-O.7 
•^•^-0.6 

o 1 4-0.8 
•^•^-0.8 

Q 1 +0.9 
-^•^-0.8 

2.6 ±0.3 

Q 2133-4311 



1-6!!;:^ 

5 7+2-0 

-3,2 

3.2)|:^ 

4 1+1-2 

4 8+2-2 

4 4+1.5 

^•^-1.4 

4.0 ± 0.7 

Q 2133-4625 






3 9+1-2 

9 4+0.7 
^•^-0.7 

9 9+1.4 
■^•^-1.4 

3.2 ±0.7 

Q 2134-4521 



2.8!];| 

Q '74-3.6 

6. /_2 9 

4 5+3-’ 

-t.j J 9 

O '74-1.4 
-1.3 

3 3+1-1 

3.6)jJ 

4.4 ±0.8 

PC 2331-1-0216 

0 . ^+0.7 

^•‘^-0.3 

-0.3 

2 7+2-1 

18.5)‘:3 

11.1)- 

3 9+1-2 

6.9)1^ 

6.1 ±0.6 

mean 

3.9 ± 0.5 

1.5 ±0.1 

1.7 ±0.2 

5.3 ± 0.9 

6.3 ± 0.8 

5.5 ± 0.6 

3.9 ± 0.5 

4.6 ± 0.6 



metallicities of the gas closely related to quasars at high red- 
shifts (z ^ 4) can be expected. To achieve the observed high 
metallicities, the single-zone and multi-zone models indicate 
evolutionary time scales for a major star formation episode of 
Tevoi - 0.5 - 0.8 Gyrs. 

Based on the overall mean metallicities of Z/Z© = 4.3 + 
0.3 of the high-redshift quasars in this study (z ^ 4) and the 
time scale for the star formation necessary to enrich the gas, the 
epoch of the first star formation can be estimated. A redshift z ^ 
4 corresponds to an age of the universe of less than ~ 1.3 Gyrs 
(Ho = 65kms“' Mpc^', Qm = 0.3, Ga = 0.7). An evolution 
time scale of Tevoi - 0.5 - 0.8 Gyrs implies that the first major 
star formation in the quasars studied here must have started at 
Z/ - 6 - 8. It is interesting to note that this is also the epoch that 
is supposed to mark the re-ionization of the universe (Haiman 
& Loeb 1998; Becker et al. 2001; Fan et al. 2002). 

Comparable high metallicities were measured for a small 
sample of quasars at z - 2.4 - 3.8 (Dietrich et al. 1999; Dietrich 
& Wilhelm-Erkens 2000), for a large sample with z - 2 - 3 
(Hamann & Ferland 1993), and for 70 quasars with z ^ 3.5 
(Dietrich et al. 2002, in prep.). In particular, there is no decline 
in metallicity from z - 2 to z ^ 4. This can be taken as an 
indication that the formation of massive spheroidal systems, 
accompanied by intense star formation, starts at z/ - 6 - 8 and 
continues until z - 2 (Madau et al. 1996; Steidel et al. 1999). 
At the end of the first major star formation episode, which lasts 
for ~ 0.5 - 0.8 Gyrs, quasar activity starts in environments that 
are already highly enriched (e.g., Granato et al. 2001; Romano 
et al. 2002). 

5. Conclusion 

We observed a sample of 11 high redshift quasars with 3.9 $ 
z $ 5 at low spectral resolution. We used several emission¬ 
line fluxes ratios involving carbon, nitrogen, oxygen, and he¬ 
lium to estimate the metallicity of the line emitting gas. To 


transform the observed line ratios into metallicities we used 
the results of detailed photoionization calculations (Hamann 
et al. 2002). The emission line ratios involving N iii] and N v 
provide generally consistent estimates of the gas metallic¬ 
ity for quasars. In particular, the results based on Niii]/Oiii], 
N v/(0 vi-hC iv), N v/O VI, and N v/C iv differ by less than ~ 
30 %. The average metallicity for the 11 high redshift quasars 
in our sample is Z/Zq - 4.3 + 0.3. We placed these re¬ 
sults in the context of chemical evolution models presented by 
Hamann & Ferland (1993) and Friaga & Terlevich (1998). For 
an evolution/enrichment time scale of approximately Temi ~ 
0.5 to 0.8 Gyrs, we estimate that the first major star forma¬ 
tion must have begun in these environments at a redshift of 
Z/ - 6 to 8, i.e., at a cosmic age of less than 1 Gyr (Hg = 65 
kms^'Mpc^', Q.M - 0.3, Da = 0.7). 

Acknowledgements. We are grateful to our colleagues J.A. Baldwin, 
G.J. Ferland, and K.T. Korista for helpful discussions. MD and FH ac¬ 
knowledge support from NASA grant NAG 5-3234 and NSF grant 
AST-99-84040 (University of Florida). JH and MD were also sup¬ 
ported by the grants SFB 328 D and SFB439 (Landessternwarte 
Heidelberg). MV gratefully acknowledges financial support from the 
Columbus Fellowship. 

References 

Allende Prieto, C., Lambert, D.L., & Asplund, M. 2001, ApJ, 556, 
L63 

Andreani, R, La Franca, F, & Cristiani, S. 1993, MNRAS, 261, L35 

Arimoto, N. & Yoshii, Y. 1987, A&A, 173, 23 

Baldwin, J.A. & Netzer, H. 1978, ApJ, 226, 1 

Becker, R.H., et al. 2001, AJ, 122, 2850 

Burstein, D. & Heiles, C. 1982, AJ, 87, 1165 

Carilli, C.L., Bertoldi, F, Omont, A., Cox, R, McMahon, R.G., & 
Isaak, K.G. 2001, AJ, 122, 1679 
Cen, R. & Ostriker, J.P. 1999, ApJ, 519, L109 
Connolly, A.J., Szalay, A.S., Dickinson, M., Subbarao, M., & Brunner, 
R.J. 1997, ApJ, 486, LI 1 



M.Dietrich et al.: High redshift quasars 


7 



redshift z 


Fig. 4. In the top panel, the average metallicity of the individ¬ 
ual high redshift quasars is shown as a function of redshift. The 
dashed line marks solar metallicity Zq. In the bottom panel the 
mean metallicity for each individual line ratio is displayed to¬ 
gether with the mean redshift range covered by the quasars for 
which the individual line ratios can be measured. 

Davidson, K. 1977, ApJ, 218, 20 

Dickinson, M. 1998, in “The Hubble Deep Field”, STScl Symp., eds. 

M. Divio, S. Fall, & P. Madau, p.219 
Dietrich, M., Appenzeller, 1., Vestergaard, M., &. Wagner, S.J. 2002a, 
ApJ, 564, 581 

Dietrich, M., Appenzeller, 1., Wagner, S.J., et al. 1999, A&A, 352, LI 
Dietrich, M., et al. 2002b, ApJ, 581, in press 
Dietrich, M. & Wilhelm-Erkens, U. 2000, A&A, 354, 17 
Fan, X., etal. 1999, AJ, 118, 1 

Fan, X., Narayanan, V.K., Strauss, M.A., White, R.L., Becker, R.H., 
Pentericci, L., & Rix, H.-W. 2002, AJ, 123, 1247 
Ferland, G.J., Baldwin, J.A., Korista, K.T., Hamann, F., Carswell, 
R.F., Phillips, M., Wilkes, B., & Williams, R.E. 1996, ApJ, 461, 
683 

Ferland, G.J., Korista, K.T., Verner, D.A., Ferguson, J.W., Kingdon, 
J.B., & Verner, E.M. 1998, PASP, 110, 761 
Friafa, A.C.S. & Terlevich, R.J. 1998, MNRAS, 298, 399 
Gaskell, C.M., Shields, G.A., & Wampler, E.J. 1981, ApJ, 249, 443 
Gebhardt, K., et al. 2000, ApJ, 543, L5 
Gnedin, N.Y. & Ostriker, J.P. 1997, ApJ, 486, 581 
Granato, G.L., Silva, L., Monaco, R, Panuzzo, R, Salucci, P., De Zotti, 
G., & Danese, L. 2001, MNRAS, 324, 757 
Grevesse, N. & Sauval, A.J. 1998, Space Sci.Rev., 85, 161 
Haiman, Z. & Loeb, A. 1998, ApJ, 503, 505 
Hamann, F. 1997, ApJS, 109, 279 


Hamann, F. & Ferland, G.J. 1992, ApJ, 381, L53 
Hamann, F. & Ferland, G.J. 1993, ApJ, 418, 11 
Hamann, F. & Ferland, G.J. 1999, ARA&A, 37, 487 
Hamann, R, Korista, K.T., Ferland, G.J., Warner, C., & Baldwin, J.A. 
2002, ApJ, 564, 592 

Hamuy, M., Walker, A.R., Suntzeff, N.B., Gigoux R, Heathcote, S.R., 
& Phillips, M.M. 1992, PASP, 104, 533 
Holweger, H. 2001, in Joint SOHO/ACE workshop “Solar and 
Galactic Composition”, ed. R.E. Wimmer-Schweingruber, AIP 
Conf.Proc. Vol.598, p.23 
Horne, K. 1986, PASP, 98, 609 

Isaak, K.G., McMahon, R.G., Hils, R.E., & Withington, S. 1994, 
MNRAS, 269, L28 

Izotov, Y.I. & Thuan, T.X. 1999, ApJ, 511, 639 
Kaspi, S., Smith, P.S., Netzer, H., Maoz, D., Jannuzi, B.T., & Giveon, 
U. 2000, ApJ, 533, 631 

Kauffmann, G. & Haehnelt, M.G. 2000, MNRAS, 311, 576 
Kormendy, J. & Richstone, D. 1995, ARA&A, 33, 581 
Kormendy, J. & Gebhardt, K. 2001, in AIP conf.proc. Vol.586, 
20th Texas Symposium on relativistic Astrophysics XIX, eds. 

J. C. Wheeler & H. Martel, p.363 
Laor, A., et al. 1994, ApJ, 420, 110 

Laor, A., Eiore, R, Elvis, M., Wilkes, B.J., & McDowell, J.C. 1997, 
ApJ, 477, 93 

Lilly, S.J., Le Eevre, O., Hammer, R, & Crampton, D. 1996, ApJ, 460, 
LI 

Madau, R, Rerguson, H.C., Dickinson, M.E., Giavalisco, M., Steidel, 
C.C., & Rruchter, A. 1996, MNRAS, 283, 1388 
Magorrian, J., Tremaine, S., & Richstone, D. 1998, AJ, 115, 2285 
Mathews, W.G. & Rerland, G.J. 1987, ApJ, 323, 456 
Matteucci, R. & Padovani, P. 1993, ApJ, 419, 485 
Merritt, D. & Rerrarese, L. 2001, ApJ, 547, 140 
Mohler, S., Seifert, W., Appenzeller, L, & Muschielok, B. 1995, in 
ESO Workshop “Calibrating and Understanding HST and ESO 
Instruments”, ed. P. Benvenuti, ESO, p.l45 
Omont, A., Cox, R, Bertoldi, R, McMahon, R.G., Carilli, C., & Isaak, 

K. G. 2001, A&A, 374, 371 

Omont, A., McMahon, R.G., Cox, R, Kreysa, E., Bergeron, J., Pajot, 
R, & Storrie-Lombardi, L.J. 1996, A&A, 315, 1 
Osmer, P.S. 1980, ApJ, 237, 666 
Padovani, P. & Matteucci, R 1993, ApJ, 416, 26 
Pagel, B.E.J. & Edmunds, M.G. 1981, ARA&A, 19, 77 
Petitjean, R, Rauch, M., & Carswell, R.E. 1994, A&A, 291, 29 
Pettini, M. 1999, in Proc.of ESO Workshop “Chemical Evolution from 
Zero to High Redshift”, ed. J. Walsh & M. Rosa, LNP, p.233 
Romano, D., Silva, L., Matteucci, R, & Danese, L. 2002, MNRAS, 
334, 444 

Salucci, R, Szuszkiewicz, E., Monaco, P, & Danese, L. 1999, 
MNRAS, 307, 637 

Savage, B.D. & Mathis, J.S. 1979, ARA&A, 17, 73 
Schwarz, H.E. & Melnick, J. 1993, The ESO Users Manual, p.24 
Shemmer, O. & Netzer, H. 2002, ApJ, 567, L19 
Shields, G.A. 1976, ApJ, 204, 330 

Steidel, C.C., Adelberger, K.L., Giavalisco, M., Dickinson, M., & 
Pettini, M. 1999, ApJ, 519, 1 

Telfer, R.C., Zheng, W., Kriss, G.A., & Davidsen, A.R 2002, ApJ, 
565, 773 

Tinsley, B.M. 1979, ApJ, 229, 1046 
Tinsley, B.M. 1980, Eundam. Cosmic Phys., 5, 287 
Tremaine, S., et al. 2002, ApJ, 574, 740 
Tresse, L. & Maddox, S.J. 1998, ApJ, 495, 691 
Uomoto, A. 1984, ApJ, 284, 497 

van Zee, L., Salzer, J.J., & Haynes, M.P. 1998, ApJ, 497, LI 






M.Dietrich et al.: High redshift quasars 


Vestergaard, M. & Wilkes, BJ. 2001, ApJS, 134, 1 
Warner, C., Hamann, R, Shields, J.C., Constantin, A., Foltz, C.B., & 
Chaffee, F.H. 2002, ApJ, 567, 68 

Wheeler, J.C., Sneden, C., & Truran, J.W. 1989, ARA&A, 27, 279 
Zheng, W., Kriss, G.A., Teller, R.C., Grimes, J.R, & Davidson, AR. 
1997, ApJ, 475, 469 



M.Dietrich et al.: High redshift quasars 


9 




4000 6000 8000 10000 


wavelength [A] 


Fig. 1. The spectra of the quasars in the observers frame. The flux is given in units of 10“'^ erg s“' A“'. Emission lines 

used in this study are labeled in the spectrum of Q 0046-293. The horizontal bars indicate the location of the continuum windows 
used to fit the continuum. 
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Fig. 3. The metallicities of the individual high redshift quasars as a function of redshift. The chemical composition which is 
derived from the emission line ratios is plotted for each ratio. 



